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RADIATIVE TRANSFER IN A 
PLANE STRATIFIED DIELECTRIC 

Thomas T. Wilheit, Jr. 


ABSTRACT 

A model is developed for calculating radiative transfer in a 
stratified dielectric. This model is used to show that the reflec- 
tivity of a stratified dielectric is primarily determined by gradi- 
ents in the real part of the refractive index over distances on the 
order of 1/10 wavelength in the medium. The effective tempera- 
ture of the medium is determined by the thermodynamic temper- 
ature profile over distances of the order 8^ where 

r 

T 

Im(n) dx = \/4 i7 
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RADIATIVE ^"RANSFER IN A 
PLANE STRATIFIED DIELECTRIC 

Thomas T. Wilheit, Jr. 


There is a burgeoning interest in microwave remote sensing of the Earth's sur- 
face. The dominant theme of both active (Radar) and passive (Radiometer) sys- 
tems is a variability of the dielectric properties of the surface as a function of 
a number of parameters of interest. It is straightforward to calculate the re- 
flectivity and microwave radiance of an interface between two homogenous di- 
electrics using the Fresnel relations (Panofsky and Phillips, 1962). However, 
quantitative understanding of many measurements requires consideration of a 
more sophisticated model than a simple planar, isothermal dielectric interface. 

It is the purpose of this document to develop a framework for consideration of 
dielectric and thermal variations in one dimension. To this end, we will develop 
a theory of radiative transfer in a planar layered medium and demonstrate that 
this layered theory is useful for consideration of continuously varying dielectrics. 

In the first section, we solve the field equations for a multiple layered dielectric 
to provide a framework for the second section wherein we solve the radiative 
transfer problem. In the third section, we calculate the thermal radiation from 
each layer and thus solve the radiative transfer problem. In the last section, we 
apply this theory to certain specific cases. 


Electric Fields within a Layered Medium 


Consider a series of N dielectric layers, each with a complex index of refraction 
Uj and a thickness 8^ . The first and last layers are semi-infinite. Electromag- 
netic waves can propagate in each layer with the form E = E^je ‘‘"j where Icj 
is the propagation vector with magnitude 277n. A , ^ is the free space wavelength, 
j specifies the layer, and r is a general position vector. E^. is an arbitrary 
electric vector with complex coefficients subject only to the constraint that 
kj • E^j = 0. We will choose the coordinate system such that the interfaces are 
normal to the x direction. In order to satisfy any (non-trivial) boundary condi- 
tions at the interfaces, we must impose Snell's law on the k vectors; 

k X X = constant 
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where x is a unit vector in the x direction. We can then further specify that 
ky - 0 with no loss of generality. We may now consider the x dependence only 
in the wave and include the y and z dependence implicity in the Eq value. 

The and - signs correspond to waves propagating in the positive and negative 
X directions, respectively. The value oi 6.^ the an^e between the direction of 
propagation and the x direction, is related among all the layers by the Snell’s 
law constraint, n^ sin - constant. 6^ is, in general, a complex number. We 
may now establish a pair of propagator functions; 

^ * (2'^i/\)J*n(x' ) cos (5(x'))dx' 

P'^(x) = e *1 


where x^ is the coordinate of the interface between the first and second layers. 

The electric fields within the j the layer propagating in either direction can be 
expressed in terms of the propagator function as this function satisfies the 
necessary wave equation at aU points. 

Defining the x origin to be at the interface of the first and second layers leads 
to a simple form for p*, the propagator at the interface between the j and j + 1**’ 
layers. 


±{2rri/\) n- 8 . 


4 J-V^Z'ii/zv^n-u-cOS 

P* = e J J 


e. i 

^Pi-i 


p 


± 

1 


= 1 


The electric fields at each interface may then be expressed in terms of this 
propagator and are shown schematically in Figure 1. It is these fields to which 
the appropriate boundary conditions must be applied in order to determine the 
electric fields. 


Any incident electric wave may be resolved into two polarizations; we will solve 
for two mutually orthogonal linear polarizations from which any arbitrary polar- 
ization may be synthesized by linear superposition. Horizontal polarization is 
defined as having the electric vector in the plane of the interfaces, whereas the 
vertically polarized wave has its electric field vector in the plane determined by 
the propagations vector and the normal to the interfaces. This definition breaks 
down at normal incidence but, by symmetry, the physical distinction between the 
two polarization breaks down as well, and the two solutions converge in the limit 
0-0. The boundary conditions to be applied are that the tangential components 
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of both the electric and magnetic fields be constant across the boundary. The 
magnetic field is related to the electric throu^ Maxwell’s equation giving: 


H = k X E/o>/x 


where co is the angular frequency of the wave and m the magnetic permitivity 
which we will assume to be the same on both sides of the boundary. Thus we 
have: 


Xx (p;e;+p“et) = Xx (pJe;,j + pTe:,j) 


and 


Xx (pjkt X e| + pTk* X E~.) = .Xx (pjk| x eJ+j + Pjk’^i x E"+j) 


Using Snell's Inw and conventional vector analysis, these may be reduced for 
horizontal polarization to: 


+ pTe: = pJe;,, + p-Ej,, 




and 


n. cos - p’EJ) = n.+j cos ^j+i(p|Et^j - PjEj+i) 

and in the vertical case to: 


and 


cos ^^(PjEj - pTE') = cos 0. + j(p.£. + j + pTE’+j) 
nj(p*E* - pTE-) = nj„(p;E‘„ - pTET,,) 


A result of our definition of the propagation function is that 

p; = (p-)-‘ 

permitting these to be further simplified 

(Pp^ E* t E- - (p*)^ E**, + e:,, 

and 

n. cos ^j[(pp2 E* - EV] = n.+j cos ^.+j[(pp2 eJ^^ - Ej.j] 
for the horizontal polarization and 


cos d. [(p|)2 Ej + E~) - cos ^j + i[(p|)^ E|+J + Ej + j] 
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and 


n. [(p|)2 E; - E-] = n, [(p*)^ - E',j] 

for the vertical polarization. 

For the N layers, we have N -1 interfaces with two boundary conditions so that 
we have a total of 2N-2 equations to solve for the 2N electric fields. We may 
also specify the incident radiation in the positive direction on the first interface 
and in the negative direction on the last interface, giving the two remaining con- 
ditions needed to determine a solution for the electric fields. Here we will 
spec; " a field of unit amplitude for the input in the positive direction and that 
no Wvive is propagating in the negative direction in the last layer. The first con- 
dition is arbitrary, but since the equations are linear, the solution may be scaled 
to whatever field is desired. The second condition requires either that there are 
no further reflection possibilities beyond the last layer considered, or that the 
intervening layers are sufficiently opaque that whatever field may exist is of no 
consequence. 

We may solve this system of equations by arbitrarily specifying E * and setting 
Ej] = 0, and then solving for E^j and E"_j using the boundary conditions for 
the last interface and continuing through each of the interfaces in turn. Since 
the resulting solution for e{ is, in general, not unity and thus does not satisfy 
the initial condition, we must exploit the linearity of the equations and appropri- 
ately scale all the E’s. 

This will enable us to solve for the net reflectively of a stratified dielectric. 

This is sulficient in the case dealing with a strong electric field so that the 
thermal radiation field may be ignored. Also up to this point, nothing new has 
been derived. This, apart from notation, is just the result obtained in chapter II 
Wait (1970). However, this notation is more convenient for the radiative transfer 
calculations which follow. 


Radiative Transfer 

Consider a black body at a thermodynamic temperature tg in the first layer 
(assumed to have n = 1 + 0 i) radiating isotropically at a rate given by Boltzman's 
law (Kraus 1966); 


B = 


2hv^ 
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where B is the power per unit time per unit cross-sectional area of radiator 
per unit bandwidth per unit solid angle, h is Planck's constant, c the speed of 
light, V the frequency, and k Boltzman's co:>stant. For frequencies less than 
about 120 GHz and terrestrial temperatures (to > 190°K),it is accurate to better 
than 0.1°K to reduce this to the RayJeigh-Jeans ^proximation; B = Zv^\st^/c^ 
so that the intensity of radiation is linearly related to the thermodynamic tem- 
perature of the radiator and may be characterized by an equivalent blackbody 
temperature, or brightness temperature; Tg - BC V2i^^k. 

Each layer absorbs a fraction k- of this radiation, where j is the index speci- 
fying the layer. If each of these layers is also at a temperature to, each must 
also radiate as much energy to the blackbody as it absorbs from the blackbody 
since they are in equilibrium. Since the directional and polarization properties 
of the black body may be varied arbitrarily, this must hold true along any propa- 
gation path and for any polarization. Since the actual radiation from any layer 
except the first is independent of the presence of a blackbody in the first layer, 
and since the amount of radiation from any layer is independent of the temperature 
of any other layer, except in so far as the temperature affects the dielectric 
properties, then the radiation propagating fn a «n\ direction, 6 , is given by 


N 

Tj(S) =2^ fi(®) ‘i +R’’(‘') 

i = 2 

where is the fraction of radiation incident on the first interface at an 

angle 0, with polarization P, which would be absorbed by the jth layer, t. is 
the thermodynamic temperature of the layer, is the reflectivity for the 

radiation (i9) incident at an angle, on the first interface. Since by con- 

servation of energy 


M 

z: 
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this reduces in the isothermal case to 

= (5) + (l-R''(c: )) to 

The crux, then, of determining the brightness tempera'^ure which would be ob- 
served from a stratified dielectric is in calculating >"he {f . } . Clearly, the value 
of f. for given layer j is (^j-i ^ i where t'j.. is tiie net electromagnetic 

energy flux entering the jth layer at the j-l’st interface and c. is the stme for 
the j + I'st layer at the j'th interface for an electromagnetic energy flux c ^ inci- 
dent on the firot interface. 
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The electromagnetic energy flow through a plane is given by the Poynting vector 
§ which is defined as; S ' ^)» ''^re denotes the real part of the 

vector product of the complex electric and magnetic fields and the asterisk (*) 
denotes complex conjugation. Since H = ic x ^/ooix^ 

s = R 

* I COfJi J 


using this relationship, it may be shown that 

/2tt = (n. cos 6 -^) I EjPj I ^ 

+ K.COj cos e.) jE./p. 

+ 2lm(n. cos ) Im(E|(E’)2 Pj /p*) 

and that 

\a.'/x8..j/27T R^(n. COS 9 .^) (E^P-.j)^ 

+ R^(n. COS 9.) |E"/p._ + 2 Im(n. cos ) Im(Ej(Ep* p..j/p*.j ) 

and since we defined the incomii^ field to have unit amplitude, 


2v COS 6*n 

?, 2 

‘ Xcofl 


so that 


S. = 

j 


1 

COS 9 ^ 


R^(n. COS 9.) lE^i^ (|Pj-iI'^ 



+ R^(n. cos 9.) i£'|^ 




,2 


+ 2 Im(n. cos 9 ^) Im 
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The two terms containing R^<n. cosd.) represent the attenuati<Hi of the forward 
and backward waves, respectively, and the term in Im(n . co80 . ) represents the 
interference between the two waves. 

In the appendix, listings of two FORTRAN IV subroutine, HPLD and VPLD, are 
presented. These subroutines solve ti^ Reid equations and calculate the thermal 
weighting functions for the horizontal and vertical polarizations, respectively. 


Sample Calculations 


As an i.pplication of this theory, consider the situation shown schematically in 
Figure 2. This represents an abrupt change in the refractive index from n| to 
n 2 , then a continuous linear change from n, to U 2 over a distance 6 and a con- 
stant index beyond that. 

This is a common situation for the air/soil interface where n^, is the refractive 
index of air and Uj is the refractive index in the deeper, homogenous part of the 
soil; nj corresponds to the surface of the soil vhich has been dried by evapora- 
tion or moistened by condensation. In Table I, the values for and n^ and corre- 
sponding moisture values for four different cases are given. These indices of 
refraction are typical of a clay loam soil at a wavelength of 19.35 GHz and were, 
along with the choice of moisture values, extracted from Schmugge et al (1974). 

It is obvious that if S is long enough (S » \).the surface reflectivity would be 
just that produced by the - n^ interface and if S« \ , the reflectivity is just 
the Uq - n^ reflectivity. The question remains of how the reflectivity behaves 
when S « \ and from what depth the thermal emission comes. 

For each of the three sets of indices, this arrangement was approximated by a 
number of layers (10, 20, 40, 80, and 160), the reflectivities and temperature 
wreighting functions wrere calculated for normal incidence and for both horizontal 
and vertical polarization at 45° incidence, and b was varied from \/lOO to 10 \. 

In each case, the result for 160 layers ivas compared with each result for fewer 
layers. The results are shown in Table 2 for the maximum difference observed 
for any S for each of the models. The reflectivity error between the result for 
80 and for 160 layers was, in no case, greater than .0017, which corresponds to 
a brightness temperature error of only about 0.5 °K, or in the case of a radar 
measurement, an error of less than 0.05 dB. The 0.5°K value is about the accu- 
racy achievable in current state-of-the-art radiometric measurements, and the 
0.05 dB value is about an order of magnitude better than the state-of-the-art in 
radar meatjurements. Thus the convergence of this stepwise approximation to 
accuracies relevant to interpretation of experiments is quite rapid. 
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In Figure 3, the reflectivity results for viewing perpendicularly to the interface 
in all four cases are plotted. In each case, for ^ > 3, the reflectivity is essen- 
tially that of the Hq - n^ interface only and for ^ A < 1/30, the reflectivity 
approaches that of the subsurface layer assy mptotic ally. There is a rather 
abrupt transition region between these values of 8 A • On the large 5 side of 
this transition, bmall local maxima and minima in reflectivity are observed. 

These correspoiid to electrical lengths for the transition region, ReCUj + mj)h/2, 
of 1/2 and 1/4 wavelength, respectivity, for the first three cases. The positions 
of the minimum and th<. maximum are reversed for th£^ fourth case. Between 
this quarter wave resonance and d « V/30, the reflectivity changes rapidly from 
a value characteristic of n ^ to a value characteristic of n 2 > We can define a 
reflectivity sampling depth for this profile, as the value of ^ at which the 
reflectivity is that \^ch woxild result from an abrupt change in the index of 
refraction from 1 to (n ^ + n2>/2. The values of 5^ are given in Table 3 for 
each the four models for perpendicular viewing and for viewing at 45® incidence 
an^e for both horizontal and vertical polarizations. The sampling depths range 
from ,073-\ to .032 The shorter sampling depths correspoiKi roughly to the 
larger indices of refraction. The sampling depths were multiplied by the average 
real part of ihe refractive index through the transition region; these results 

are also given in Table 2. These values show little variation from a minimum of 
0.132 to a maximum of 0.162 for the first three cases. They are somewhat 
smaller for the wet surface case. 

It Is of particular interest to note the lack of significant angular dependence. 

This occurs because the direction of propagation through the medium remains 
nearly normal to the interface even for oblique incidence as a result of Snell’s 
law. The decrease of the angular effect with increasi’ig dielectric constant is 
consistent with this interpretation. 

It is also possible to calculate the average depth at which the thermal radiation 
upwelling from the soil originates. We will call this the thermal sampling depth, 

, and it is given by 


S 


T 



rf 
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where Xj is the depth of the ith layer and is the weighting function for that 
layer as previously defined. It is readily shown for a uniform dielectric 
(n, = Dj) that 


^ 4-rr Im(n) 
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This may be generalized as an approximation to to aid in understanding the 
physics of the problexu. We define S' such that 


f*T 

4n I Im(n) dx = A. 

•'0 

hi Fig^e 4f the values of (solid line) and (dashed Tine) are shown for tlm 
four perpendicular viewing cases. In cases two and three* the dashed line coin- 
sides with the solid line and thus is not visible. Note that for even the wettest 
cases, the thermal sampling depth is greater than the reflectivity sampling depth 
and for the drier cases is more than an order of magnitude bigger than . 

The close agreement between and Indicates that the thermal sampling 
depth is primarily a function of integrals over the imaginary part of the index 
of refraction. The minor descrepancies v^ch do appear result primarily from 
ignoring absorption of the reflecte wave within the medium. These calculations 
were repeated with all the imaginary parts of the indices of refraction reduced 
by a factor of ten. The thermal sampling ecpths were drastically different, as 
would be expected, but the reflectivity sampling depth changed by less than 1!5% 
in all cases. This indicates that the reflectivity is determined, for the most 
part, by gradients in tue real part of the refractive index. 

Thus, for understanding radiometric measurements of surfaces, two sampling 
depths must be considered. The reflectivity charactei .zes changes in the rcial 
part of the refractive index over a depth of about l/lO to 1/7 wavelengths in the 
medium. The thermal radiation is characterized by a different depth which 
is generally larger and is substantially determined bv the imaginary part of the 
index o* refraction. 
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Appendix I 

FORTRAN IV SUBROUTINES 

The two subroutines presented here, HPLD and VPLD, calculate the electric 
fields and thermal weighting functions for the Horizontally and Vertically polar- 
ized layered dielectric, respectively. The number of layers to be treated is N, 
the angle of incidence in radians, THETA, and the wavelength is XLAM. It re- 
quires 4 complex and 1 real single precision vectors of length N. The complex 
numbers are signified by $. $N is the vector of complex indices of refraction 
of the layers. $SEP and $EM are computed by the subroutine and are the elec- 
tric field solutions for the positive and negative directions, respectively. DEL 
is the vector of layer thicknesses; the first and last of these are not iised. $P 
is a work vector. It contains tie values of Pj during computations. The values 
of the thermal wei^ting function are rytuiaed in the real part. The power re- 
flectivity is considered to be the weighting fimction for the first layer. 
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SUBROUT FNE HPLU« tN, »EP . tEM.OEL* SPt Nf T hET A t XL AM | 

IMPLICIT complex I»» 

DIMENSIONSNI 10), SEP ( 10),SP( 10 ) , SEM ( 10) « DEL ( 10 ) 

S=SINITMETA) 

$pm = c l,, 0 .) 

NLsN-1 

on 10 I=?,NL 
NMAX=I 

SS=SN( I )»S/SN( I ) 

SC=CSORTH 1 . ,0. )-SS*SS) 

ARG=DEL( I )*?*3.1A159?/XLAM 
SARR=?.*ARG#SN< I )*SC#(0.,1.) 

SP( I)=CEXP( SARG)*SP( I-l ) 

IFICABSCSPI !) ).LT. 0.0001 ) GO TO 20 
10 CONTINUE 
PQ CONTINUE 

SEP(NMAX)=( l.,0.) 

SFM(NMAX)=(0.,0.) 
on AO JJ=2,NMAX 
J=NMAX-JJ+1 
SSJ=SNI I )*S/SN( J) 

SCJ = rSORT( ( I .,0. )-SSj=»*2 ) 

SSJP1 = SN( 1)<'S/SN( J+1) 

SCJP1=CS0RT( ( I . ,0. )-SSJPl*#2) 

SA = 2*SN( J )*SCJ/( SN( +SN(J+1 )*SCJP1 ) 

SR=( SN( J )*$CJ-SN(J+1 )«SCJP1 )/( (SNIJ |#SCJ + SN( J+ 1 ) *$C JP 1 ) *SP ( J ) ) 
SFP(J)=SEP( J+1.)/SA+SB#SEM( J+1 )/SA 
SFM(J)=SEM( J+I )+( SEPIJ+1 )-SEP( J) )*SP( J) 

40 CONTINUE 
fX=SEPI 1 ) 
on 45 J=1,NMAX 
$FPU)=SEPC J)/SX 
SEMI J)=SEM( J)/SX 
45 CONTINUE 

no 50 J=NMAX,N 
50 SPI J ) = (0. ,1 . )/() ,E50) 

LL=NMAX-1 
DO 60 .IJ = 1,LL 
J=NMAX-JJ +1 
SS=SIN(THETA)/SN(.I) 

SC=CSORT( ( 1. ,0. )-.5S<^«2 ) 

R = CARS( SPI J ) ) 

5=CARSI SPIJ-1 ) ) 

F2=I S-R I’^CAKSI SEPUIl'i'O^+n./K-l . /S ) »»CABS ( SEM I J) )=>=»2 
np=e?*RFALI SNIJ ) *SC ) /COS! THETA ) 
sxp=sEPi J )«cnNjG( SEMI .1 n 

X^2.«A IMAGI SN'IJ )#SC/CnSITHETA ) ) # I A IMAG I SX P#SP IJ- 1 ) /C Afl S I SP IJ- I ) ) ) 
C-AIMAGI SXP^SPI J )/CABS ( SPIJ) ) ) ) 
f)P = nP-X 

SPI J )=CMPLXIOP ,0. ) 

60 CONTI N'.;f 

R = CAHSI '•EH' ! > -’i-i^wEAl. I SNI I ) ) 

SPI 1 ) =r.MP‘.x : p } 

RFTUHN 

FNO 


ORIGINAL PAGE IS 
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SIIRKHUT 1NE VPLn(«'N,%E>'»».EMtnFL»S>M,N.THFTAtXLftNi) 

implicit complex cm 

DIMENSIONtNC 1 0 ). *E P ( 1 0 )♦ 1.P (]<»)» S>EM ( lO )» Ot L <1 O ) 

S=S1N( TMETA ) 

I. P( I ) = ( l.,0. ) 

Nl. = N-l 

DO IM I=i>,NL 
NMAX= F 

$S = tM( 1 )#S/ J>N( I ) 

SC = CS(.)«T( ( 1. tO. > 

APr, = nEL( n*?-3. 1^1‘5 4X/XLAM 
SA«r,= ?.*A«C*tN( ! )*tc*(n. ♦ 1 . ) 

J. PF I )=CFXP(1.AP(, )s;%P( I-] ) 

FEFCAHSF tP( I ) F.LT.O.nOOl ) CP If' ?0 

in CONTI MUF 
?n CONTFNUF 

*EP( MMAX ) = ( 1 , ,0. ) 

SFMF (viMAX ) = (0. »0. ) 

Dn C»n J.| = ?,NMAX 

J = NMAX-.I J+1 

SSJ = ^N I 1 )*S/ tMF J ) 

SC J = CSO«T ( ( 1 . ,0. )-S‘SJ=f*? > 

SSJP1 = SM( 1 ) VS/SN.M .1 + 1 ) 

SCJHl=CS!.)kT( (l.,0.)-SSJPl--?F 
sf)=:^.-sN( J )=? sr J 
SA = SN( J)s=SCJPl + SN( .)+! )*SC.J 
t,R = $NF ,) )^sr,JPl-SN( J+1 )-SCJ 

SFP( .) ) =SA=?SEP( J+1 ) /SO +SH-.SEM( J+1 )/(M)-SP( J) I 
SR = SM{ .1+1 )/SN« J ) 

SFM{.I)=SKX^SFM(J+1 ) + FSEF(J)-SFP(.l + l F-SKl-SPFJ) 

40 CONTINOF 
SX=SFP ( 1 ) 
nn 4S .1 = 1 ,NMAX 
SFP( .J ) =SE P( J » /SX 
SFM( ,I)=SFM( .) )/%X 
4R CONTINUE 

no SO ,I = NMAX.M 
bO ,SP( ,) 5 = ( n, t 1 . ) / ( 1 .E‘iO ) 
i L= npax-1 
no +,0 ,l.l = l,LL 
,) = NMAX- l.l +l 
SS = S I K ( theta ) /%^.( .1 ) 

-lOCAC KT ( ( I . )-SS--P ) 

k=CAHS( SP( ,1) ) 

S = C.A‘<S( SP( J-1 ) ) 

t-p = (S--<>vrARS(SFP(.J) )v=.^;? + ( 1. /P-1. /SI -C/.K<s ( SFF ( J ) F -'P 
np = E F A|. { +.N( .) ) - SC ) / COS ( 1 hF T A I 
^,vp=hFP( .) )-CnfJ.IM SE .'( .1 ) ) 

X = ;> . V A I |'A(.( SN( ,1 ) -S( /COS ( 1 »<F [ A ) ) V ( A ! I'i A(. ( S XP-tP ( .) - 1 ) / cars ( S p ( .) - 1 ) ) I 

C- A I M AC. ( S X P - S P ( .) ) / (, Af« S ( S P ( .1 ) ) ) ) 

IIP^ >p - X 

Sp ( .1 ) = f PL X ( lip , n. ) 

RO (, liA'IINUf 

P = C A P S ( SI K ( 1 ) ) •V - P P A l_ ( S N ( 1 ) ) 

f ( 1 ) = I V I X ( P , i ) . ) 
k h Ti ■ 
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Figure 1. Schematic representation of reflected and 
transmitted waves at the interface between the j'th 
and j + l^et layer. 
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TRANSITION ZONE THICKNESS 8 

Figure 3. Power reflection at normal incidence as a function of transition 
zone thickness. The circled numbers by the curves give the model number 
(see Table 1). 
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Table 1 

Indices of Refraction and Moisture Values 


Model 


% Moisture 
Dry Weight Basis 

nj 

% Moisture 
Dry Weight Basis 

n 

1.8 + O.OSi 

3 

2.2 + 0.25i 

12 


2.4 + 0.28i 

15 

3.1 + 0.601 

20 


3.8 + 0.901 

22 

4.0 + 0.98i 

25 

n 

4.0 + 0.98i 

25 

1.8 + O.OSi 

3 


Table 2 

Convergence of Multilayer Approximation 


Model 

Viewing 

Angle 

Polarization 

Reflectivity Error 

®^AYERS 

40 

LAYERS 

^ Layers 

LAYERS 

1 

0° 

- 

.0003 

.0010 

1 

RR 

1 

45° 

H 

.0005 

.0012 

1 ! 

mmm 

1 

45° 

V 

.0002 

.0005 

^R 1 1 

RB 

2 

0° 

- 

.0005 

.0012 


.0116 

2 

45° 

H 

.0005 

.0017 


.0111 

2 

45° 

V 

.0003 

.0010 

^R K 

.0065 

3 

0° 

- 

.0001 

.0003 

^R m ^R 

.0017 

3 

45° 

H 

.0001 

.0002 

^R m ^R 

.0017 

3 

45° 

V 

.0001 

.0003 

^R m ^R 

.0016 

■■ 

0° 

- 

.0017 

.0051 


.0245 


45° 

H 

.0017 

.0052 


.0297 

■1 

45° 

V 

.0015 

.0047 


.0264 
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Table 3 

Reflectivity Sampling Depths 


Model 

Viewing 

An^e 

Polerization 


<“r> 

<^,>8, A 

1 


— 

.066 

2.0 

.132 

1 

45° 

H 

.071 

2.0 

.142 

1 

45° 

V 

.073 

2.0 

.146 

2 

0° 

- 

.056 

2.75 

.154 

2 

45° 

H 

.056 

2.75 

.154 

2 

45° 

V 

.059 

2.75 

.162 

3 

0° 

- 

.040 

3.90 

.156 

3 

45° 

H 

.040 

3.90 

.156 

3 

cn 

o 

V 

.040 

3.90 

.156 

4 

0° 

- 

.033 

2.90 

.096 

4 

45° 

H 

.033 

2.90 

.096 

4 

45° 

V 

.032 

2.90 

.093 
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